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We use a recently developed self-consistent GW approximation to present systematic ab initio 
calculations of the conduction band spin splitting in III-V and II-V zincblende semiconductors. The 
spin orbit interaction is taken into account as a perturbation to the scalar relativistic hamiltonian. 
These are the first calculations of conduction band spin splittings based on a quasiparticle approach; 
and because the self-consistent GW scheme accurately reproduces the relevant band parameters, 
it is expected to be a reliable predictor of spin splittings. The results are compared to the few 
available experimental data and a previous calculation based on a model one-particle potential. We 
also briefly address the widely used k-p parameterization in the context of these results. 
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In this letter, we apply a recently developed all- 
electron quasiparticle (QP) self-consistent GW method 
(QPscGW) 0, S 13 to determine the spin splitting of 
conduction bands in III-V and II- VI zincblende semicon- 
ductors as a consequence of the spin-orbit coupling. This 
quantity is emerging as a parameter of great importance 
in spintronics applications, in part because it determines 
spin lifetimes, but also because it may be possible to 
exploit the splitting to induce spin current without an 
external magnetic field through the Rashba effect. 

Information about this important quantity is very 
sparse. On the theoretical side, precise QP levels for 
the occupied and first few unoccupied states are required. 
We show here that QPscGM^ accurately reproduces mea- 
sured fundamental and higher-lying gaps in all the semi- 
conductors studied, as well as other key parameters (ef- 
fective mass, d band position, etc.), and is therefore a 
suitable vehicle for accurate determination of the split- 
ting. 

Owing to a lack of inversion symmetry, the conduction 
band of a zincblende semiconductor is split by spin-orbit 
coupling. Dresselhaus Q showed that near the minimum 
point Pg, the splitting can be described by the opera- 
tor Hi) = • Bcff, where a are the Pauli matrices and 
Big = 2jki{kf^-^ - kf_^_2), i = {x,y, z}. Hb is written in 
this form to suggest that its effect is equivalent to the ac- 
tion of a k-dependent effective magnetic field Bcff, which 
is anisotropic in k and is proportional to fc'^. 7 is a con- 
stant that depends on the bulk properties of the material. 
Provided the Larmor precession frequency under i?cfj(k) 
is less than the inverse of momentum relaxation time, 
small random changes in k will cause the polarization of 
an injected current to gradually diminish. This is the 
spin scattering mechanism of Dyakonov and Perel (DP) 
I3; and it is widely accepted that for a wide range of 
temperatures and moderate concentration of impurities, 
DP is the dominant spin relaxation mechanism in the 
conduction band of III-V and II- VI semiconductors. The 
DP spin scattering rate can be estimated theoretically 
once the value of 7 is known . 



In low dimensional structures the splitting scales in 
proportion to fc, owing to confinement effects, which 
causes the DP mechanism to increase significantly. More- 
over asymmetry in the leads or structure of a device re- 
sults in an additional spin splitting effect, proposed first 
by Rashba and Bychkov This contribution from this 
"structural inversion asymmetry" depends on the bound- 
ary conditions: it can interfere constructively or destruc- 
tively with the Dresselhaus term depending on device 
eometry and k vector. Based on this, a few authors 
1,0 proposed that growth of quantum wells in certain 
crystallographic directions can result in long spin coher- 
ence times, if the ratio of Dresselhaus to Rashba terms 
is close to unity. Measurements of these two contribu- 
tions suggest that they can be comparable 0, ^3 ■ Also 
giant spin relaxation anisotropy can occur in strained 
bulk zincblende semiconductors owing to destructive in- 
terference between the Dresselhaus and strain-induced 
spin splittings 0- 

To evaluate 7 theoretically, it is essential to use an ab 
initio approach rather than a parameterized method such 
as the k-p theory because the latter includes many exter- 
nal parameters whose values can only be reliably deter- 
mined through fitting to an ab initio approach. However, 
the ab initio approach based on the usual local density 
approximation (LDA) is no better. This is because 7 
depends strongly on key features of the band structure, 
such as effective mass, that the LDA cannot reliably pre- 
dict. As is well known, it is difficult to choose an ap- 
propriate empirical adjustment to the LDA potential so 
that the energy band structure of the lowest few conduc- 
tion bands is reliable. This is particularly true in more 
complex cases such as quantum wells where experimental 
data is not available. 

QPscGW is a method to determine a (nonlocal but 
static) hermitian one-particle hamiltonian Hq in a self- 
consistent way With the usual GW approximation, 
we can calculate the self-energy S(w) as a functional of 
the starting trial Hq. Here S(ijj) can be expanded in the 
basis of the eigenfunctions {V-'kn(r)} of Hq, where k is 
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TABLE I; Important band parameters for the III-V semiconductors. Eg and E'q are the energies of the first two conduction 
bands at the F-point, and are defined in the text; Eg is the fundamental gap when it differs from Eo- /m is the conduction 
band effective mass at F. Energies are in eV; 7 is in eV-A^. Quantities in rectangular brackets were computed from Eq. 2 with a 
chosen to reproduce the experimental fundamental gap i?g. Quantities in parenthesis and curly brackets correspond respectively 
to experimental data, and calculated data from Ref.ll2L Spin orbit splittings Aso and Ago do not change significantly whether 
the QPscGW, the scaled QPscGW or the LDA potential is used. 
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the wave vector and n is the band index. Then we define 
the static self-energy as 

-D / / I ^(fkn) + 5](ekn') I ; \ /i\ 
Sknn' = Re^V'knl ^ [V'kri') (1) 

where, ekn denote eigenvalues, j,nd Re means to take the 
hermitian part. By using this Skrm' instead of the usual 
LDA exchange correlation potential V^J-'^, we can iterate 
to self-consistency, i.e. the Sknn' generated by Hq is 
identical to the I]k„„' that enters into Hq. We explicitly 
include core contributions to E. 

With the addition of local orbitals, we include in the 
MTO basis all atomic-like states that are within ~ ±2 Ry 
of the Fermi energy: thus both the 3d and Ad levels are 
included for GaAs, and so on. As QPscGVF gives the 
self-consistent Hq at the scalar relativistic level, we add 
the spin-orbit operator Hso — L ■ S/2c? to Hq as a per- 
turbation (it is not included in the self-consistency cycle) . 
Thus the band structure is calculated from a one-body 
Hamiltonian Hq + i?so- We also considered a 'scaled E' 



with 

Ho. + Fso = ^LDA + (1 - - V^^^) + Hso. (2) 

(Note that i/LDA and V^^^ are also determined from 
the self-consistent density given by Hq). Tables I and 
II present QPscGVF results (a = 0) for III-V and II- VI 
compounds. We also present results in a few cases where 
a is chosen to reproduce the experimental fundamental 
gap Eg at OK, which we do to address the strong sensi- 
tivity of 7 to Eg. The error in Eg from QPscGM^ the- 
ory, while small, is comparable to Eg itself for small-gap 
systems InAs and InSb, resulting in large errors in the 
effective mass m^. Evidently once (S — V^^^) is scaled 
so Eg is reproduced, falls in good agreement with the 
experimental mass. The magnitude of 7 depends on the 
hybridization between the lower conduction band and the 
p-like conduction and valence bands. The hybridization 
depends on the position of the first two conduction band 
levels Eq = E{T%)-E{Tl) and E'q = E{r^) - E{T^s). (In 
the absence of spin orbit coupling, Fg has Fi symmetry; 
FojFS and FS,FS are the valence and conduction levels 
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with symmetry.) and is reflected in the magnitude of 
effective mass m^. Also important is the strength of spin- 
orbit interaction at anion and cation sites, reflected in the 
spin orbit splittings Ago and AgQ of Fj'g and FJg . With 
the exceptions of AgQ in InP and Eq in AlP, all param- 
eters are in good agreement with available experimental 
data. The experimental AgQ seems low in InP and incon- 
sistent with the trends described below. The discrepancy 
in £^0 for AlP is also curious as it falls well outside the 
error usually expected from QPscG'W^:2j; also its better 
known fundamental gap agrees well with QPscGM^. It 
is likely that Eq was misidentified in the (30 year old) 
experiment. 

As the valence band has mostly anion character and 
the conduction band mostly cation character, one would 
expect that Ago is controlled mainly by the anion, while 
AgQ is controlled mainly by the cation. Comparing com- 
pounds with the same cation but different anions, or vice- 
versa, this trend becomes evident. For example, for GaP, 
GaAs and GaSb A^q is 0.158, 0.174, and 0.196 eV re- 
spectively; while on the other side for AlAs, GaAs and 
InAs A^Q is 0.031, 0.174 and 0.429 eV. This trend is 
found in all materials presented in Tables I and II. 

Fig. n] shows the conduction band spin splitting ver- 
sus wave number for the entire [110] line. For sufficiently 
small k it varies as as expected from the k-p analysis; 
7 was obtained by fitting a polynomial to the calculated 
splitting at k points near F. The sign of 7 was deter- 
mined from the eigenvectors and is defined according to 
the most accepted convention {AE ~ i?(S4) — £(£3)). 

The splitting is rather complex, arising from several 
factors such as how the F^g and F^g states admix with 
the conduction state as k moves away from F. Still, in- 
spection of Table 1 and Fig. 1 reveals some trends: 
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(i) 7 scales approximately in proportion to 1/Eq, at 
least in the narrow-gap case, as can be seen by compar- 
ing scaled {Ha + Hso) and unsealed {Hq + i?so) results 
for InAs and InSb. This is expected by the k-p theory. 
The scaling only slightly perturbs the potential; the band 
structure stays approximately constant except for small 
shifts in the lowest conduction band that substantially 
affect the gap and the conduction band shape near the 
F-point (and the valence band that couples to this band). 

(ii) 7 tends to increase as the spin-orbit coupling param- 



eter A =3 ((||)v|^^so|(||)c) becomes more negative. 



FIG. 1: Conduction band spin splitting in [110] direction. 
GaAs, GaSb, InP, InAs and InSb are plotted for the scaled 



A becomes more negative with increasing anion mass, 
while the reverse is true when cation mass increases. 

When anion mass increases, these tendencies add con- 
structively and 7 increases monotonically. When cation 
mass increases, these tendencies interfere with each other, 
and the behavior is more complex. 

Table I compares our calculations of 7 with the work 
of Cardona et al Those authors added empirical 

pseudo-Darwin potential shifts to make their band struc- 
ture match the observed bandgap. These shifts result 
in a reasonable and qualitatively correct band structure, 
though band parameters have varying degrees of accu- 
racy (e.g. deviated from experiment as much as 
35%[l3)- The splittings we obtain are in qualitative, 
and often quantitative agreement with Cardona's calcu- 
lations; e.g. AE changes sign at approximately the same 
fc-points in GaP and GaAs. 

What limited experimental data that exists is also 
shown in Table I. For GaAs [H, GaSb H^, InP [13, 
7 1 was inferred from an estimation of the spin relax- 
ation rate ascribed to the DP spin relaxation mechanism. 
This method does not resolve the sign of 7. For InSb 7 
was also extracted from electric-dipole-magnetic-dipole 
interference at far infrared frequencies |20j . which is a 
more direct measure of the splitting and resolves the sign. 
This is probably the most accurate measurement of 7 
in any III-V semiconductor. The measurements of 7 in 
Refs. [ll[l3 are rather indirect and are based on extract- 
ing simultaneously and accurately spin and momentum 
relaxation times. Also there can be some uncertainty 
related to the type of momentum scattering mechanism. 

Some care must be taken to determine 7 in InSb. 
Owing to the small gap, the splitting be gins to devi- 
ate from a k^ dependence at very small k [21|. For the 
electron densities used in that work (7Ve=1.6x lO^** to 
4.7 X 10^^ cm~^) the Fermi wave number kp ranged be- 
tween 0.9 X 10~^ and 2.7 x 10~^ a.u.~^ At this k, we 
find a small deviation from the k^ dependence ('^6%). 
However at larger doping the renormalization is signifi- 
cant: at = 10 X 10""^ a.u.~^ the splitting is reduced by 
40% (i.e. AEhk^ = 0.6) in InSb, by 20% in InAs, but 
is negligible in the remaining semiconductors. 

In Fig. 121 we show the LDA and QPscGW^ spin split- 
ting together with the corresponding bands for GaAs. 
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TABLE II: Important band parameters for the II- VI semi- 
conductors. Symbols are the same as in Table I. Experimental 
data are taken from Ref. fTJ unless noted otherwise. 



ZnS ZnSe ZnTe CdTe 

'E'o Im 2^94 2^45 1.84 

(expt) (3.78) (2.82) (2.394) (1.606) 

E'q-Eq 4.59 4.52 2.91 3.57 

Aso 0.084 0.408 0.884 0.829 

(expt) (0.065)° (0.4) (0.91) (0.81)° 

Ago 0.140 0.133 0.110 0.305 

A" -HO. 16 -0.15 -0.44 -0.28 

'mJi 0.186 0.131 0.113 0.100 

(expt) (0.184) (0.13) (0.13) (0.096)'' 

7 -0.48 -fl.29 -H13.3 +8.5 

{Ref. 12} {+1.6} {+11-7} 



"From Ref. na| 
''From Ref. 121 



Roughly speaking, the spin splitting of the lowest con- 
duction band is induced by hybridization with the cation- 
like conduction band and the anion-like valence bands of 
p symmetry. These two bands contribute to 7 with op- 
posite sign. Dispersion in the relevant bands on the [110] 
line differs significantly in the LDA from the QPscGM^ 
case; but the relative difference is largest near T point 
where the LDA overestimates the sp hybridization. Con- 
sequently the relative difference in LDA and QPscGW^ 
conduction band spin splitting is largest near F: 7^^^ is 
+121 eV-A'^ in GaAs, which is very different from the 6.4 
or 8.5 eV-A'^ in TablcHJ Near X where the splitting is lin- 
ear in k, LDA and QPscGVF splittings are very similar. 
This trend is common to all materials tested. 

Until now, several authors have relied on the k-p ap- 




FIG. 2: (a) The LDA (red dashed line) and QPscGW (solid 
blue line) upper valence and lower conduction bands, (b) 
The LDA (red dashed line) and QPscGW (solid blue line) 
conduction bands spin splitting in [110] direction. 



proximation to estimate the value of 7 [l^ I21I . Gen- 
erally speaking, the "downfolded" k-p Hamiltonian (i.e. 
with the far band contributions folded in) can contain 
all possible independent parameters consistent with the 
crystal symmetry at the P point. However, it is not possi- 
ble to determine these parameters just from experiments. 
Thus it is necessary to assume a number of constraints so 
as to reduce the number of the parameters. For example, 
a usual approximation is that the Luttinger's "residual 7 
parameters" for the unoccupied p bands are taken to be 
zero I23I . However, constraints vary from paper to paper 
[1^ l23l ]. and the spin splitting one obtains depends on 
the choice of constraint. In principle, all parameters can 
be calculated from the GW hamiltonian without mak- 
ing any additional assumptions, in which case the down- 
folded Hamiltonian should reproduce all the GW QP lev- 
els exactly near P. However, the parametrization of k-p 
theory is based on the Luttinger-Kohn expansion of the 
V'kn |2J], where the periodic part of '0kn at k ^ is ex- 
panded by Tpon- This is neither convenient nor suitable 
for the parametrization of ah initio eigenfunctions. It is 
on the other hand rather straightforward to expand all 
the required matrix elements in the basis of the ab initio 
ipkn that respect the crystal symmetry. 

In conclusion, we employed a recently developed 
QPscGW^ to make ab initio predictions of the conduc- 
tion band spin splitting in III-V and II- VI semiconduc- 
tors. (In a few small empirical adjustment to the 
self-energy was also needed.) Such an approach is essen- 
tial to reliably determine the splitting; moreover it builds 
a framework that will enable builders of model hamilto- 
nians to reliably incorporate the effect without the need 
for additional assumptions. 
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